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The Effect of 53 /im IR Radiation on 18 cm OH Megamaser 

Emission 
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ABSTRACT 

OH megamasers (OHMs) emit primarily in the main lines at 1667 and 1665 
MHz, and differ from their Galactic counterparts due to their immense luminosi- 
ties, large linewidths and 1667/1665 MHz flux ratios, which are always greater 
than one. We flnd that these maser properties result from strong 53 /im radia- 
tive pumping combined with line overlap effects caused by turbulent linewidths 
~ 20 kms~^; pumping calculations that do not include line overlap are unre- 
liable. A minimum dust temperature of ~ 45 K is needed for inversion, and 
maximum maser efficiency occurs for dust temperatures ~ 80 - 140 K. We find 
that warmer dust can support inversion at lower IR luminosities, in agreement 
with observations. Our results are in good agreement with a clumpy model of 
OHMs, with clouds sizes ^ 1 pc and OH column densities ~ 5xl0^^cm~^, that is 
able to explain both the diffuse and compact emission observed for OHMs. We 
suggest that all OH main line masers may be pumped by far-IR radiation, with 
the major differences between OHMs and Galactic OH masers caused by differ- 
ences in linewidth produced by line overlap. Small Galactic maser linewidths 
tend to produce stronger 1665 MHz emission. The large OHM linewidths lead 
to inverted ground state transitions having approximately the same excitation 
temperature, producing 1667/1665 MHz flux ratios greater than one and weak 
satellite line emission. Finally, the small observed ratio of pumping radiation to 
dense molecular gas, as traced by HON and HCO+, is a possible reason for the 
lack of OH megamaser emission in NGC 6240. 
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Introduction 



OH megamasers are extremely luminous extragalactic OH maser sources that have 
isotropic luminosities a million or more times greater than their Galactic counterparts. OHMs 
are found only in the nuclear regions of luminous (LIRGs) and ultraluminous (ULIRGs) in- 
frared galaxies, where intense star formation is occurring. These masers emit primarily in 
the main lines at 1667 and 1665 MHz, although weaker emissi on in the satellite lines at 1612 
and 1720 MHz has also been detected (IBaan fc Haschicklll987l ). In addition to their immense 
luminosities, OHMs differ from their Galactic counterparts in their very large linewidths and 
their 1667/1665 MHz flux ratios. Galactic OH mainline masers in star forming regions tend 
to have linewidths < 1 kms~^ and 1667/1665 MHz flux ratios less than one. In contrast, 
OHMs have overall linewidths > 100 kms~^ and their 1667/1665 MH z flux ratios are always 
greate r than one. Although this ratio varies, the typical ratio is ~ 5 (ILonsdald |2002| ) . iBaan 
(119891 ) found a clear separation between OH emitters and absorbers based on the IR proper- 
ties of the host galaxies. The OHMs have larger IR luminosity (>Cir) and tend to be warmer. 
OHMs may occur fo r lower £tr as long as the dust is warm enough. The Arecibo OHM 
survey conducted by iDarling fc Giovanellil (120021 ) targeted LIRGs and in c reased the number 
of OHMs to ~ 100. That survey substantiated the conclusions of iBaanI (119891 ) and showed 
that the fraction of OHMs in LIRGs is an increasing function of far-IR (FIR) luminosity 
and color, with most of the detections found in warm ULIRGs. These observations strongly 
suggest that the masers are pumped by the intense FIR radiation in the OH pump lines at 
53 fim and 35 fim. Suppor t for FIR purnping is provided by observations using the Infrared 
Space Observatory (ISO). ISkinner et al.l (119971 ) observed the 35 /im pumping transition in 
Arp 220 and concluded that absorption in this line alone could power its maser emission. 
Signi ficantly, OH and water megamasers have not been observed in the same galaxy (ILo 
2005! ) even though both are often found in the same Galactic star forming regions. 



The standard OHM model was introduced by iBaanI (119851 ). who proposed that the 
maser emission is produced by low gain unsaturated amplification of background radio con- 
tinuum. This proposal was confirmed by the comprehensive study of OHMs performed by 
Henkel fc Wilson] (jl990l ). They found good agreement with observations if the maser tran- 
sitions have approximately equal excitation temperatures, giving a 1667/1665 MHz optical 
depth ratio of ~ 1.8. However, subsequent VLBI observations revealed compact maser emis- 
sion o n parsec scales with a mplification factors > 800 and very large 1667/1665 MHz line 
ratios ( Lonsdale et al.l 1998 ). A surprise of these observations was tha t linewidths rem ained 
large (tens of kms~^) even on the smallest observed angular scales (ILonsdald |2002| ) . The 
VLBI observations led to the suggestion that there are two different modes of maser oper- 
ation — low gain, unsaturated amplification responsible for the diffuse maser emission and 
high gain, saturated emission producing the compact sources. It was also suggested that the 
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two classes may have different pumping mechanisms, with the diffuse emission p umped by 
IR ra diation and the compact masers by a combination of colhsions and radiation (ILonsdale 
20021 ). However, an extension of the standard model is able to explain b oth the diffuse an d 
compact emission from IIIZw35 by assuming a clumpy maser medium (jParra et al.l l2005l ) . 
Each cloud generates low gain, unsaturated emission, and strong compact emission occurs 
when the line of sight intersects more than one cloud. A similar model has been used to 
explain the megamaser emission from IRAS 17208-0014 (IMomjian et al.ll2006l ). The clumpy 
maser model also explains the observation that comp act masers are always found embedded 
in the diffuse emission and do not occur in isolation (ILonsdald |2002| ) . 



An observational difficulty is that the nearest OHMs are located at distances of ~ 100 
Mpc and thus cannot be probed nearly as accurately as Galactic masers. The study of 
OHMs thus relies heavily on analyzing their global properties. It is important to first ex- 
plain the global properties of OHMs and then construct more comprehensive models for 
those sources which have been observed in more detail u sing VLBI techniques. Previous 
OHM models (IHenkel fc WilsonI Il990l : iRandell et al.l Il995l ) were performed before the dis- 
covery of the compact maser sources. Although the clumpy maser model seems to provide 
a phenomenological explanation to the compact and diffuse maser emissions, it is not yet 
backed by pumping calculations. The time is right to develop a pumping model that can 
explain both the compact and diffuse emission using physical conditions that are consistent 
with those expected in the maser regions. This is the aim of the present paper. In section 
2 we summarize what is known concerning the physical conditions existing in the maser re- 
gion. In section 3 we explain our model and present the results of our calculations. Section 
4 contains a summary and discussion. 



2. Physical Conditions in the Maser Region 



Calculating maser intensities requires knowledge of the physical conditions existing in 
the maser region. The important parameters are the FIR radiation intensity, OH column 
density, gas density and temperature, and linewidth. The physical conditions in the nuclear 
regions of OHMs are significantly different from those found in normal galaxies due to the ex- 
treme amount of mas sive star formation. Although hidden from view at visible wavelengths, 
Lonsdale et al.l (120061 ) have observed numerous radio point sources in the nuclear region of 
Arp 220 and conclude that they are radio supernovae. They deduce a supernova rate of ~ 4 
per year, yielding a per-volume rate that is orders of magnitude larger than Galactic. These 
radio sources have recently been o bserved at other r adio wavelengths and the SN nature of 
the sources has been strengthened (jParra et al.ll2007l ). That study also found that an unusu- 
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ally large fraction of these supernovae are highly luminous, suggesting that stellar formation 
and evolution may be much different in these regions compared to normal galaxies. The 
large amount of massive star formation results in the intense FIR radiation and turbulence 
present in the region of maser production. 



2.1. FIR Radiation Intensity 



The intense FIR radiation field present in the nuclei of LIRGs and ULIRGs is the 
dominant pump for the OHMs. Estimates of the effective dust temperature and opacity 
have been calculated in numerous studies. The IR AS derived col o r temperatures of 74 
OHMs range from 40 - 90 K with a median of 60 K jKlockneil l2o"o2l ) . IYuu fc Cariuil J2002h 
find best fit dust temperatures in the range 59 - 72 K for the OHMs Arp 220, Mrk 273 and 
Mrk 231. These are global es t imate s and may be influenced by cooler dust lying outside 
the maser region. ISoifer et al.l (119991 ) performed high spatial resolution observations of the 
nuclear region of Arp 220 from 3 to 25 fim and conclude that the region is complicated with a 
range of source temperatures and sizes. They find that emission at wavelengths longer than 
20 fim corresponds to a dust temperature of ~ 85 K. This result is consistent with a detailed 
theoretical analysis of radiation pressure-supported starbursts, which finds the effective dust 
temp erature to be ~ 90 K for the optically thick starburst disks of ULIRGs ([Thompson et al. 
20051 ). A possible upper limit to the dust temperature of ~ 100 K ma y be provided by th e 



lack of observation of the 6.7 GHz methanol maser emission in OHMs (jParling et al.l 120031 ). 
This methanol transition is often found to be a strong maser in association with OH masers 
in Galactic star forming regions. Modeling of the 6.7 GHz maser fin ds that it is radiatively 
pumped and requires dust temperatures > 100 K (ICragg et al.ll2005l ). 



The dust opacity in OHMs is not constrained as tightly by observations , and the esti- 
mated opacity depends on the wavelengths being examined. iLisenfeld et al.l (120001 ) discuss 
a number of arguments that suggest visual opacities of ord er 1000. However , other studies 
estimate much smaller visual opacities in the range 20-200 (ISoifer et al.lll999l : iHinz &: Rieke 
20061 ). It should be noted that the values of dust temperature and opacity that best fit 
the observed spectral energy distributio ns are not indepe ndent and there is considerable 
degeneracy between the two parameters (IBlain et al.ll2003l ). 
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2.2. Gas Temperature and Molecular Density 



Gas temperature and density are needed to calculate collision rates. Although OHM 
pumping is dominated by radiation, collisions need to be included in the calculation of 
the OH level populations. It has also been sugges ted that the co mpact masers may be 
pumped by a combination of collisions and radiation (jLonsdaldl2002l ) . The temperature and 
density can be es t imate d using high density tracer molecules such as HCN, CS and HCO^. 
Cernicharo et al.l (120061 ) observed transitions in HCN and HNC in Arp 22 and their analysis 



predicts densities in the range 10^-10^ cm ^ and gas temperature ~ 40 K. lGreve et al.l (120061 ) 



observed Arp 220 using numerous molecular lines including those of HCN, CS and HCO"*". 
They deduce densities in the range 10^-10^ cm~^ and gas temperatures ranging from 50 K 
to 70 K. Further constraints on these parameters are provided by the lack of detection of the 
22 GHz water maser co mbined with the detection of the 183 GHz water maser in Arp 220 
(jCernicharo et al.ll2006l ). The 22 GHz maser is coUisionally pumped and requires relatively 
high densities and temperatures for its production. The 183 GHz maser is produced at 
much lower densities and temperature. ICernicharo et al.l (l2006l ) perform calculations that 
show that strong 183 GHz emission without 22 GHz emission requires densities < 10^ cm~^ 
and gas temperatures < 100 K. 



2.3. Linewidth 

Linewidth usually affects the level populations only as part of the combination tt-oh-R/ 
the column density per velocity interval; neither rioH? R nor AV enter separately. However, 
overlap between FIR rotational lines of OH introduces the linewidth as an independent 
parameter with additional, separate effects on the level populations. Line overlap affects 
photon trapping since a photon emitted in one transition may be absorbed in a different 
one. In addition, line overlap can significantly affect the absorption of the pumping radia- 
tion in the overlapping transitions. The lowest 24 hyperfine levels of OH are shown in Figure 
[1] The level structure is more complicated than is found for the typical diatomic molecule. 
There are two rotational ladders, and each rotational level is split into 4 sublevels due to the 
effects of lambda doubling and hyperfine splitting. The hyperfine separations differ in the 2 
halves of the lambda doublets leading to line separations as small as .1 kms~^ between FIR 
rotational lines. Figure [1] shows that the asymmetry is especially large for the transitions 
involving the ^Hi/2 ladder. 

An important feature of OHMs is their very large linewidth. Single dish observations of 
OHMs reveal linewidths of hundreds of kms~^. Much of this linewidth is due to large scale 
motions of individual maser clouds. However, VLBI imaging finds that linewidths are still 
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tens of kms~^ even for parsec size regions (iLonsdald |2002| ) . These observations indicate the 
presence of a l arge amount of superson ic turbulence and are in agreement with the theoretical 
predictions of iThoinpson et al.l (120051 ) . who estimate turbulent linewidths of ~ 50 kms~^. 



Downes &: SolomonI (119981 ) performed a detailed analysis of the central region of ULIRGs 



using CO interferometry and found that the contribution of local turbulence t o the overall 



linewidth is > 30 km s ^ . Detailed modeling of the megamasers in IIIZw35 by iParra et al. 



(120051 ) found best agreement with observations using a cloud internal velocity dispersion of 
20 kms~^. These large linewidths lead to a significant amount of overlap in the OH FIR 
transitions. 



3. Modeling and Results 



Although all four ground state transitions have been observed as masers in OHMs, the 
satellite lines at 1612 MHz and 1720 MHz are much weaker than the main lines at 1667 MHz 
and 1665 MHz. Satellite line maser pumping follows naturally from the OH level structure, 
but main line pu mping require s an asymmetry in the excitation of the two halves of the 
lambda doublets (lElitzurl 119921 ). Possible asymmetric pumps include radiative pumping by 
warm dust or overlap between the OH FIR rotational lines. In the Appendix we discuss FIR 
pumping in detail and show that line overlap completely dominates all other suggested FIR 
pumping mechanisms of main line masers. 

We solve for the populations of the 24 levels shown in Figure [H These include all the 
levels with direct radiative connection to the ground state, where the masers are produced. 
We have checked that increasing the number of levels makes no difference on the outcome. 
The stati s tical equilibrium equations are set in the escape probabihty approach using the 
Capriottil (119651 ) escape probabihty. This is equivalent to a slab approximated with a constant 
source function, and should provide reasonable estimates of the overall emission from the 
region. Further discussi on of this technique and compar ison with exact radiative transfer 
solutions can be found in lElitzur fc Asensio Ramod (120061 ). As noted in §2.3^ the OH energy 
level structure necessita tes inclusion of the effects of line overlap. We use the method of 
Lockett fc Elitzurl (119891 ) to treat the local overlap of the FIR rotational lines. This procedure 
is based on the slab escape probability method and accounts for overlap effects by numerically 
integrating the amount of absorpti on in the overlapp i ng lin es assuming Gaussian line profiles. 
This method is similar to that of iGuilloteau et al.l (jl98ll ). except they assume a spherical 
geometry as opposed to a slab. 



We envision the maser source embedded in a dusty environment that provides the pump- 
ing IR radiation. Assuming the dust temperature T^i to be approximately uniform across the 
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maser cloud, the local radiation field can be described with 



/, = 5,(T,)(l-e-^-) 



1^ 



Here ByiTd) is the Planck function and Ty is the frequency dependent dust opacity across 
the region where its temperature can be described by the single value T^. Note that Ti, is 
the dust optical depth over the entire dusty region and can be larger than it s value across 
themaser cloud. The spectral shape of t^^ is found using the cross sections of iDraine fc Lee 



Although radiation is the dominant pumping mechanism, collisions need to be i n clude d 
in the calculations and we utilize the most recent collision rates from lOffer et al.l (119941 ) . 
These rates are considerably different for collisior is with ortho and pa ra H21 and the pumping 
could depend sensitively on the ortho-para ratio (ILockett et al.lll999l ). This ratio is uncertain 
in molecular clouds. H2 is believed to form with the ortho-para ratio ~ 3 and the ratio 
approaches th ermal equilibrium at a rate depending on the chemical history of the cloud 
(iFlowerl Il990l ). We performed calculations ranging from pure ortho to pure para and the 
results did not vary significantly. Our nominal choice for the ortho-para ratio is one, which 
is the equilibrium value for T ~ 77 K. 

We explored a wide range of the physical parameters and found a large volume of phase 
space where the main lines are inverted. Table [T]lists the values of the parameters we consider 
"standard" . We now present our results, varying the model parameters one at a time from 
these nominal values. 
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Dust temperature, Td 60 K 

Dust optical depth at visual 100 

Linewidth (FWHM), 20 kms-^ 

OH column density, noH-R 6x10^^ cm~^ 

Gas temperature 50 K 

Gas density 10'^ cm~'^ 



Table 1: "Standard" model parameters chos en to be representative of individual clouds in 
the clumpy OHM model of (jParra et al.ll2005l ). Figures show the effects of varying individual 
parameters from these nominal values, as indicated in the corresponding captions. 
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3.1. Effects of dust temperature and opacity 



Figure [2] plots ttie contours of 1667 MHz negative optical depth versus dust temperature 
and opacity. A minimum dust temperature of about 45 K is needed for maser production, 
and maximum maser inversion occurs for dust temperatures between 80 - 140 K for the 
expected range of dust opacity. When the dust opacity becomes very large, the dust radiation 
approaches that of a blackbody and the maser begins to thermalize. Beyond a visual opacity 
of about 300 there is little dependence of maser strength on dust temperature. There is 
also a minimum FIR pumping flux needed for maser production, that is a function of dust 
temperature. The pumping flux in equation 1 is determined by the product of two factors: 
the Planck function, which depends on dust temperature, and the self-absorption factor, 
which depends on the dust opacity at the pumping frequency. The FIR flux is increased 
by increasing either the dust temperature or opacity. Figure 3 shows the minimum 53 /xm 
flux needed for inversion as a function of dust temperature. A larger flux is needed for 
a lower dust temperature, indicating that the FIR spectral shape is important for maser 
pumping. Thus a smaller Ci^i can support i nversi o n, as long as the temperature i s warm 
enough. This agrees with the observations of iBaaru (119891 ) and lDarling fc Giovanellil (120021 ). 
Our analysis also shows that a minimum dust temperature is needed for inversion. Below 50 
K, the Planck function at 53 /im decreases very rapidly as dust temperature decreases. The 
flux can be increased by increasing dust opacity, but it cannot exceed blackbody emission 
and the needed flux can not be produced for dust temperatures below about 45 K. 

The pumping is primarily by the 53 /im lines, with the 35 /im transitions in a secondary 
role. Increasing the number of OH energy levels in the model calculations, main-line inversion 
was not produced until the ^Hi/2(^ = 3/2) levels, at 53 /xm, were included. The 53 /xm lines 
dominate the pumping, but the maser strength increased further when the four levels of 
^Hi/2(J = 5/2), at 35 /xm, were added. We flnd the opacity of the 53 fim lines t o be about 
5 times those of the 35 /im lines. This agrees with the study of iHe fc ChenI (120041 ) who 
compared the line strengths in Arp 220 and found the amount of absorption in the 53 /im 
line to be about 6 times that of the 35 /im line. Results were hardly affected by inclusion of 
the ^H3/2(J = 7/2) levels, as expected since these levels are radiatively disconnected from 
the ground state. 

An important observed correlation is the increase of £oh with Ci^. Early OHM ob- 
servations revealed an approximately quadratic relation between the two (IBaaru Il989l ). but 
subsequent studies using a much larger sampl e and taking into account Malmquist bias es- 
tablished the flatter dependence Lqh oc L};| (jParling fc Giovanellil 120021 ). A linear relation 
is expected because low gain ampliflcation of background radio continuum implies that Cqh 
is proportional to the radio luminosity, which in turn is proportional to £ir over the relevant 
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luminosity range (jParling fc Giovanellill2002l ). This linear proportionality will not be altered 
by the expected variations among different sources in OH covering factor of the radio con- 
tinuum, wh ich wil l only induce a scatter around it. Since the observed correlation is steeper 
than linear, iBaaru (119891 ) conjectured that maser optical depth Tmaser also varies linearly with 
the pump luminosity £ir,. However, our calculations show that, depending on the different 
parameters, Xmaser in fact can ei ther increase or decrease with the IR luminosity. Indeed, in 
their extensive study of OHMs, iHenkel fc WilsonI (119901 ) did not find a correlation between 
£iR and Tmaser- A possiblc explanation for the increase of £oh in excess of l inear proportion 
to Cm may be provided by the study of HON emission from LIRGs by iGao fc Solomon 
( I2OO4J ). They find that Ci^ oc Cucn and conclude that the amount of dense (rij^j > 10^ 
cm~^) molecular gas, which is traced by the HCN, is linearly proportional to £ir. It is thus 
reasonable to assume that, together with the increase of dense molecular gas, the number of 
OH maser clouds is also increasing with £11^. 



Ld (120051 ) pointed out that, based on its FIR color and Cm the ULIRG NGC 6240 



should be an OHM (see also figure 2 of lBaanlll989l ). A possible reason why it is not is its 
low ratio of C m to dense molecular gas. T his gas is traced by HCN and HCO^ molecules, 
and Table 1 in iGracia-Carpio et al.l (120061 ) shows that the ratios Ar/'^hcn and Ar/'Chco+ 
are both less than half in NGC 6240 than in every OHM, seven in all, included in their 
sample. Thus the FIR pump rate per molecule in NGC 6240 is much lower than in the 
typical OHM. Therefore, a minimum ratio of pump i ng rad iation to dense molecular gas may 
also be a prerequisite for producing OHM. iDarlingI (120071 ) has recently conducted a detailed 
study that shows a clear separation between OHMs and OH absorbers based on the amount 
of dense molecular gas. 



3.2. Effects of collisions 

We found no conditions where collisions alone produce the observed maser emission. 
When combined with a strong FIR radiation field, collisions tend to weaken and then ther- 
malize the masers. Increasing the collision rate by increasing either gas temperature or 
density tends to reduce the maser strength. An increase in the collision rate also increases 
the dust temperature required for inversion — increasing the density from 10^ to 10^ cm~^ 
raises the minimal dust temperature from about 45 K to 50 K. These results argue against 
collisional pumping as a source of the compact masers. Another argument against collisional 
pumping is the absence of the 22 GHz water maser combined with the presence of the 183 
GHz water maser in Arp 220. We calculated the level populations for ortho and para water 
for conditions representative of the OHM region, and our results are consistent with those of 
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Cernicharo et al.l (120061 ) . The presence of the 183 GHz maser combined with the absence of 
the 22 GHz maser imphes densities < 10® cm~^ and temperature < 100 K. These conditions 
are not conducive to colhsional pumping of the OHM. 



3.3. Effects of Line Overlap 

We find that the effects of hne overlap dominate the main line inversion even at linewidth 
as small as 0.5 kms~^. In the Appendix we discuss FIR pumping in detail and show that 
line overlap completely dominates all other suggested FIR pumping mechanisms of main line 
masers. Calculations that do not include overlap effects cannot produce reliable results for 
FIR pumping. 

Figure H] shows the effects of line overlap on the 1667 MHz maser strength and the 
1667/1665 MHz optical depth ratio. Because optical depth is determined by the column 
density per velocity interval, hqiiR/ AV remains a scaling variable, as is evident from the 
figure's right panels; they show that the pumping results display similar behavior for the 
different linewidths when plotted in terms of this variable. Line overlap effects add AV as 
another independent parameter of the pump scheme, with striking effects on the inversion: 
the model results for the ratio R = r(1667)/r(1665) split into two groups, with a sharp 
transition between the two. Starting from AV = 2 kms~^, a decrease of 1 kms~^ has little 
effect on R while an increase of 1 kms~^ increases R significantly, especially at the lower 
end of column densities. Linewidths > 10 kms~^ have R ~ 1.8 independent of AV. All 
ground state hyperfine transitions are then inverted with approximate l y the same excitation 



temperature, an equality that was also noted by iHenkel &: Wilson! (jl990l ). In contrast, 
linewidths < 2 kms~^ produce R < 1, again displaying behavior that is nearly independent 
of AV. Further discussion of this behavior is provided in the Appendix. 

Our results can explain why the flux ratios observed for OHMs differ radically from 
those found for their Galactic counterparts. Galactic interstellar masers tend to have small 
linewidths and the 1665 MHz line is usually the strong est, with the strength of the 1667 



MHz line often comparable to that of the satellite lines (IGaume fc Mutell 119871 ). However, 
there are exceptions where the 1667 MHz line or one of the satellite lines is the strongest. 
Our results show that this is the expected behavior for AV < 2 kms~^. In contrast, in 
OHMs the 1667 MHz line is always the strongest and the satellite lines are much weaker by 
factors > 10, the expected behavior when AV > 10 kms^^. Significantly, main line masers 
in OH/IR stars tend to have larger linewidths than Galactic interstellar masers and their 
1667/1665 MHz flux ratios tend to be larger than one. 
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We find that all ground state transitions are inverted with approximately equal excita- 
tion temperatures when AV > 10 kms~^. This is con sistent with the OHMs global fluxes 
determined from observations (IHenkel fc Wilson! llQQOl ). The optical depth of the 1667 MHz 
line is then ~ 1.8 times that of the 1665 MHz line and ~ 9 times that of the satellite lines 
and the fluxes for the OH maser amplified background continuum become functions of op- 
tical depth alone. Figure [5] displays the amplified 1665, 1667 and 1720 MHz line fluxes. 
Because of its stronger amplification, at a given instrumental sensitivity the 1667 line will 
be detected over the largest portion of a source. The 1667/1665 and 1667/1720 flux ratios 
increase as the 1667 MHz maser gets stronger. Small optical depths give a 1667/1665 MHz 
flux ratio of ~ 2 and the typically observed ratio of 5 implies the magnitude of r(1667) is 
~ 3.5. Because the 1720 MHz is so weak compared to the 1667 MHz line, its detection is 
limited. Figure E] plots the track of 1667/1720 MHz flux ratio versus the 1667/1665 MHz 
flux ratio a s r(1667) is varied . Also show n are t he d ata for four relative ly nearby OHMs, 
taken from iBaan et al.l (119891 ). iBaan et al.l (119921 ) and iMartin et al.l (119891 ). The agreement 
is surprisingly good considering the large uncertainties in the data. 



3.4. Clumpy OHM 



The global properties of OHMs are consistent with the standard model of low gain am- 
plification of background radiation with the maser transitions having approximately equal 
excitation temperatures. Compact, bright maser emission is explained by the clumpy na- 
ture of the amplifying medium. Each maser cloud produces low gain, unsaturated emission. 
Compact emission would be observed when the line of sight intersects a number of maser 
clouds. This r nodel has been used to give very good agreement with the maser emission 
from HIZw35. iPihlstrom et al.l (120011 ) observed HIZw35 using intermediate and high reso- 
lution, and detected both diffuse and compact emission. They proposed a model based on 
a clumpy maser med ium located in a rotating ring. This model was extensively improved 
by lParra et al.l (120051 ) . who used Monte Carlo simulations to account for the overlap of indi- 
vidual maser clouds. The clumpy maser model gives results that are in excellent agreement 
with the observed spectra and with interferometric observations on all angular scales. The 
compact maser emission arises at the tangents to the ring where cloud overlap is greatest 
in both position and velocity. The model of HIZw35 finds best agreement with observations 
using maser clouds having a size < .7 pc, magnitude of r(1667) ~ 1.5 and velocity dispersion 
of ~ 20 kms~^. A similar, but less detailed, model h as since been develop ed to explain the 
compact and diffuse emission from IRAS17208-0014 (IMomjian et al.ll2006l ). 



Our "standard" parameters in Tabled] were selected to match the HIZw35 model results 
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of iParra et al.l (120051 ) . Figure [7] is an expanded view of the relevant portion in figure [Hand 



shows that the magnitude of r(1667) is about 1.5 at an OH column density of ~ 5x10^^ cm~^. 
With the model density of 10^ cm~^, this implies an OH abundance of ~ 2xl0~® for the 
cloud size of 0.7 pc. Calculated OH column densit ies are typicall y in the range 10^^ to 10^^ 



cm ^ for both PDRs (IHartquist fc Sternberglll991l ) and C-shocks (IWardlelll999l ). Larger OH 



column densities may be obtained through other mechanis ms such as the photod issociation 
of water produced from the evaporation from grain mantles (IHartquist et al.lll995l ). However, 



even that mechanism has an upper limit of ~10^^ cm~^. Figure H] shows that very large 1667 
MHz optical depths can be produced, but only if AV is < 20 kms^^ or the OH column 
exceeds ~10^^ cm^^. Neither condition is likely to be met in OHMs. Maser amplification 
by single clouds is weak because linewidths are large (> 20 kms~^) and OH columns are 



limited to 5, 10^^ cm ^. 



The clumpy maser model can also be used to explain the lack of compact maser emission 
from Mrk 231. This i s the only OHM examined with VLBI not to exhibit compact emission 
(ILonsdale et al.ll2003l ). Furthermore, the 1667/1665 M Hz flux ratio has been found to be 



2 on angular scales ranging over 3 orders of magnitude (IRichards et al.ll2005l : lKlockner et al. 



20031 ). These observations are consistent with the main lines having similar excitation tem- 
peratures while undergoing low gain am plification. Speculativ e explanations for the lack of 
compact emission have been proposed (ILonsdale et al.l l2003l ) . but a simpler e xplan ation 



is based on th e geom etry of a rotating clumpy maser disk. IRichards et al.l (120051 ) and 
Klockner et al.l (|2003l ) each propose that the masers are produced in a rotating disk that 
is more inclined to the line of sight than is found for OHMs such as IIIZw35. The length 
along the line of sight through the disk is thus smaller than if t he disk were observed edge 
on, resulting in little cloud overlap and no compact emission. IRichards et al.l (120051 ) esti- 
mate maser cloud sizes ~1 pc with the m agnitude of r('1667 ) about 1. These estimates 
are comparable to the values suggested by iParra et al.l (120051 ) in their detailed analysis of 
niZw35. 



4. Summary and Discussion 



OHMs are characterized by their extreme luminosity, 1667/1665 MHz flux ratio greater 
than one, and large linewidths. Our detailed pumping analysis of OHMs shows that the 
above properties are a natural outcome of an intense FIR pump combined with line overlap 
effects due to large turbulent linewidths. The 53 /xm lines provide the primary pump, and a 
minimum dust temperat ure of ~_45 K is n eeded for main line maser pr oduction. In agreement 
with the observations of lBaanI (119891 ) and lDarling Sz Giovanellil (120021 ). we find that a smaller 
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£iR can support inversion, as long as the temperature is warm enough. We find no conditions 
where colhsions can produce the observed maser emission, and when combined with a strong 
FIR radiation field, collisions tend to weaken and then thermalize the masers. These results 
are consistent with the detection of the 183 GHz water maser combined with the lack of 



detection of the 22 GHz water maser in Arp 220 (jCernicharo et al.ll2006l ). The 22 GHz 
maser is collisionally pumped and requires relatively high densities and temperatures for its 
production, while the 183 GHz maser is produced at much lower densities and temperature. 
The large OHM linewidths produce significant overlap among the the 53 fim pump lines 
resulting in nearly equal negative excitation te mperatures for all grou nd state transitions, in 
agreement with observations (figure [6l see also lHenkel fc Wilsonlll990l ). The 1667/1665 MHz 
fiux ratios are thus greater than one and only weak satellite line emission is produced. Our 
pumping results are in agreement with those required by clumpy maser models that explain 
both the diffuse and compact maser emission from OHMs such as IIIZw35. Finally, the small 
observed ratios of £ir/£hcn and >Cir/>Chco+ relative to those found in OHMs is a possible 
reason for the lack of OH megamaser emission in the ULIRG NGC6240. A minimum ratio 
of pumping radiation to dense molecular gas may be a prerequisite for becoming a OHM. 

Excited state transitions at 5 cm and 6 cm have also been observed in OHMs. However, 



these lines have only been detected in absorption. iHenkel et al.l (119871 ) observed absorption 
in the three 6 cm ^Hi/2(^ = 1/2) transitions in five OHMs, a nd absorption has been de- 
tected also in one 5 cm '^Il^/2{J = 5/2) transition in Arp 220 (IHenkel et al.lll986l ). These 



observations are surprising since masing in the 5 cm and 6 cm lines has often been observe d 
in conjunction with Galactic ground state OH maser sources (ICesaroni fc Walmsleylll99ll ). 
Our pumping models of ground state maser emission in OHMs usually predict absorption in 
the 5 cm lines. However, the 6 cm lines are weakly inverted with opacities a hundred times 
smaller than the 1667 MHz line. Our calculations find that absorption in the 6 cm lines 
requires smaller dust temperatures than exist in the maser regions. We expect the absorp- 
tion to occur outside the maser regions where the dust temperature is lower. Understanding 
the absorption at 5 cm and 6 cm is an important problem that we plan to treat in a future 
paper. 



4.1. OH main-line maser pump 

Our results suggest that all OH main-line masers could be pumped by the same mech- 
anism: far-IR radiation. The major differences between Galactic star forming regions on 
one hand and evolved stars and OHMs on the other can be attributed to differences in line 
overlap effects. As shown in §3.3[ linewidth is the most important factor in determining the 
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1667/1665 MHz flux ratio and thus is the dominant factor causing the extreme differences 
between the flux ratios of the two classes. In Galactic star forming regions the observed 
linewidths are < 1 kms~^ and the 1665 MHz maser is usually strongest, with the 1667 MHz 
line often comparable in strength to the satellite lines. In OHMs the linewidths are large (> 
10 kms~^) and the 1667 MHz is always stronger than the 1665 MHz line, with the satellite 
lines much weaker. Our results show that such large linewidths produce approximately equal 
negative excitation temperatures for the four ground state transitions, resulting in the 1667 
MHz line being stronger than the 1665 MHz line and much weaker satellite lines. Smaller 
linewidths tend to produce a stronger inversion in the 1665 MHz line. In evolved stars, 
the main-line maser velocities are larger than those found in Galactic star forming regions 
and they also tend to have 1667/1665 MHz flux ratios > 1. Thus the linewidth could be 
the main reason why OH main-line masers display differences in spite of sharing a common 
pump mechanism. 

Partial support by NSF grant AST-0507421 (M.E.) is gratefully acknowledged. 

A. Main-line Pumping by FIR Radiation 

Pumping of the ground state main line masers at 1667 and 1665 MHz requires asym- 
metry in the excitation of the two halves of the lambda doublets. The main ingredients of 
the pump analysis are as follows: 

1. Most molecules are in the ground state, from which they are excited into higher levels 
and cascade back to the ground. 

2. Downward transitions are due to radiative decays. For each level on the ^Hi/2 ladder, 
the strongest decay is always down the ladder. Thus the cascade proceeds primarily 
along the ^Hi/2 ladder, until the flnal decay '^Ui/2{J = 1/2) — » ^H3/2(J = 3/2). Pumps 
that rely on cross-ladder transitions are thus inherently weak. 

3. Line overlap modifies the fundamental radiative transition rates, making the FIR pump 
cycle more conducive to asymmetries. 

The last point is especially important. In the absence of overlap, the escape probability 
approximation for the radiative rate is Af3, where A and f3 are, respectively, the transition 
A-coefficicnt and the escape probability. For optically thick transitions, Aj3 A/t, and the 
radiative decay rates then become independent of line strength. The distinctions between 
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different transitions disappear in this case, malfing it difficult to maintain any asymmetry 
required for inversion. This is no longer true in the presence of line overlap. The transition 
rate now depends on the strengths of all overlapping lines in a complex way that involves 
also the degree of overlap. The distinctions between different transitions are maintained over 
a much larger range of column densities, creating the possibility for stronger maser effect. 

FIR main line pumping is dominated by line overlap effects, primarily in the 53 fim 
pump. On the ^113/2 ladder, the hyperfine splits are roughly equal in the two halves of the 
A-doublet. But on the ^Ili/2 ladder, the hyperfine split of the upper half of each A-doublet 
is larger than in the lower half by factors of 4-6 (see figured]). These uneven splits create an 
asymmetry in the overlap of the pump lines connecting the ground state to the rotational 
states of the ^Ili/2 ladder. In contrast, there is essentially no overlap asymmetry in the 
image transitions for the 120 fim pump lines connecting the ground state to the the first 
excited state in the ^113/2 ladder. The 53 fim pump lines have special properties because 
they populate the ^Ili/2{J = 3/2) state, whose quantum numbers directly mirror the ground 
state. This places it in a special position during both the excitation and cascade phases. 
First, the 53 /xm pump cycle is the only one in which each upper level connects to the 
two ground state levels in the appropriate A-doublet; in all other cases, one of the excited 
rotational levels has only one transition to the ground because of the dipole selection rules. 
This doubles the number of potential overlaps compared to the other cycles, creating more 
asymmetry during the excitation phase of the pump cycle. Second, barring cross-ladder 
transitions, this is the only pump cycle that shuffles the molecules between the two halves of 
the ground-state A-doublet. In contrast, both the 80 fim and the 35 fim pump cycles return 
the molecules into the very same A-doublet from which they started. 

The large efficiency of the 53 /im pump is produced by asymmetric overlap effects on both 
the absorption rates out of the ground state levels and the retu rning decays. The inverting 



effects of line overlap on the abs orption rates were first noted by iBurdyuzha fc Varshalovich 



(I1973I ) and applied to OHMs by IBurdyuzha &: Vikulovl (Il990l ). They found that the smaller 



amount of overlap between the 53 /xm pump lines out of the lower A-doublet of the ground 
state reduces the shielding from the pumping dust radiation a nd thus increases the pumping 



out of the lower A-doublet relative to that from the upper. IBurdyuzha fc Vikulovl (Il990l ) 
used a fixed Doppler width of .7 kms~^ and found the 1667/1665 MHz flux ratio to be ~ 1, 
in agreement with our results. 

The asymmetries of the radiative excitations disappear when AV > 10 kms~^ because 
there is then complete overlap of the 53 //m lines. However, the inversion is sustained at 
larger linewidths due to overlap effects in the cascade back to the ground. Overlaps among 
the ^Hi/2(J = 3/2) ^Hi/2(J = 1/2) cascades are particularly important since they are the 
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primary decay routes out of the upper rotational state of the 53 /im pump. In addition, the 
overlap asymmetry between the 2 halves of these A-doublets is extreme. The separations 
between the transitions connecting the lower doublets range from .5 to 2.5 kms~^, while the 
corresponding transitions between the upper doublets have separations ranging from 3.3 to 
15 kms~^, a factor of 6 larger. The line overlap asymmetry in these transitions remains 
important for very large linewidths. The overlap of these lines significantly enhances the 53 
fim pump efficiency and is necessary to invert the ground state A-doublets for linewidths > 
10 kms~^. 

When the linewidths exceed ~ 40 kms~^ the asymmetries in the cascade transitions 
di sappear, too, and the g r ound state lines thermalize. It may be noted that the calculations 



of iBurdyuzha &: Vikulovl (Il990l ) did not produce inversions at AV^> 5 kms^^ because they 
neglected overlap effects in the decays. In contrast, we include all possible overlaps, which 
essentially involve every FIR transition when AV > 10 kms^^. 
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Fig. 1. — The 24 hyperfine levels of OH used in our calculations. A-doublet and hyperfine 
splitting are not to scale. Note the very uneven spacings in the lambda doublets of the ^Ili/2 
ladder, which cause extremely asymmetric line overlap between the FIR lines involving these 
states. 
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Fig. 2. — Contour plots of the inverted 1667 MHz optical depth versus dust temperature 
and optical depth. All other model parameters are as listed in Table [T] 
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Fig. 3. — Minimal 53 /im intensity required for inversion — cooler dust requires larger 
minimal flux. At a given dust temperature, the plotted intensity corresponds to the minimal 
Ty that produces maser inversion (see eq. 1). Results are normalized to the minimal intensity 
when = 120 K (1.08 xlO~^^ ergcm"^ Hz~^ ster~^). Model parameters are in Table [H 
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Fig. 4. — 1667 MHz optical depth and 1667/1665 MHz optical depth ratio for different 
linewidths. Left panel is plotted versus OH column density. Right panel is plotted versus 
OH column density per velocity interval. All other model parameters from Table [TJ Note 
that linewidths < 2 kms~^ lead to 1667/1665 MHz optical depth ratios < 1, while linewidths 
> 3 kms~^ cause this ratio to increase dramatically, becoming ~1.75 for linewidths > 10 
kms~^. 
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Fig. 5. — Variation of flux with 1667 MHz optical depth for the three indicated masers 
when all ground-state transitions have the same excitation temperature (the situation when 
Av > 10 kms~^; see §3.3p . The 1667/1665 and 1667/1720 flux ratios increase as the 1667 
MHz maser gets stronger. The typically observed 1667/1665 MHz flux ratio of 5 implies 
r(1667) ~ —3.5 and the 1720 MHz line will always be weak compared to the 1667 MHz line. 
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Fig. 6. — "Color-color" diagram for OH maser lines. Filled squares show data points for 
four OHMs. The solid line is the track of flux ratios of the amplified background continuum 
when all lines have the same excitation temperature. The 1667 MHz optical depth varies 
along the track from —0.5 at the lower left to —4 at the upper right. 
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Fig. 7. — 1667 MHz maser optical depth for an OHM cloud using the model parameters 
from Table [H This is a zoom in on the relevant portion of figure 4. The optical depth of 
—1.5 for OH column density of ~ 5xl0^^cm~^ is in good agreement with the requirements 
of the clumpy model of OHM emission in IIIZw35 (see §3.4p . 



